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NUMERICAL  SOLUTION  OF  HYPERSONIC  FLOW  NEAR  LEADING  EDGE  OF  A 
FLAT  PLATE 

Fu  Dexun 

Beijing  Institute  of  Aerodynamics 
Abstract 

This  paper  uses  a  one  step  difference  scheme  with  second- 
order  accuracy  in  space  for  the  numerical  solution  of  Navier- 
Stokes  equations  as  brought  forth  in  Reference  [ 1 ]  to  numer¬ 
ically  calculate  the  transonic  and  hypersonic  interaction 
problems  near  the  leading  edge  of  a  flat  plate.  The  scheme  is 
discussed  and  analyzed  with  simple  model  equations.  We  used 
this  method  to  calculate  the  interaction  problems  of  the  leading 
edge  of  a  flat  plate  with  M<o=20,  Tw/To=0.06,  Re^  L=4xl03  and 

M-,  =3,  T  =T  .  and  Re_=103. 
f  W  0  ®  L 

The  calculation  results  accord  well  with  the  test  results  and 
with  these  results  we  can  more  clearly  understand  the  flow 
characteristics  of  the  interaction  near  the  leading  edge  of  a 
flat  plate. 

I.  Introduction 

What  we  are  interested  in  here  is  a  steady-state  solution 
for  the  interaction  problems  near  the  leading  edge  of  a  flat 
plate.  Starting  from  the  use  of  non-steady-state  partially 
differential  equations  to  solve  its  steady-state  problems  is 
the  simplest  iteration  method.  Being  that  the  steady-state 
solution  is  necessary,  the  solution  process  then  does  not  re¬ 
quire  a  solution  corresponding  to  the  original  non-steady-state 
partially  differential  equations  but  it  is  only  necessary  that 
the  stabilized  solution  approach  the  solution  of  the  steady- 
state  problems.  This  gives  the  structural  difference  scheme  a 
certain  amount  of  flexibility.  This  paper  proposes  an  improved 
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iteration  scheme  based  on  the  non-steady-state  equations.  It 
is  a  one-step  scheme  and  does  not  approximate  the  original  non¬ 
steady-state  differential  equations.  Yet  its  stabilized  solution 
has  second  order  accuracy  in  space.  It  can  be  known  through 
analysis  that  the  stable  conditions  of  this  scheme  are  more  re¬ 
laxed  than  those  of  most  explicit  schemes  and  it  is  even 
simpler  than  the  commonly  used  two-step  difference  scheme.  This 
paper  uses  this  method  to  calculate  the  interaction  problems  of 
the  leading  edges  of  dual  flat  plates. 

In  the  last  several  years,  many  theories  have  appeared  and 
a  great  deal  of  research  has  been  done  on  the  hypersonic  flow 
around  the  leading  edge  of  a  flat  plate  [1,8].  Test  results 
show  that  in  the  hypersonic  flow  around  a  pointed  flat  plate, 
there  exists  a  very  small  non-continuous  flow  region  near  the 
leading  edge  of  the  flat  plate.  This  region  uses  the  quantity 
of  the  molecular  average  free  path  and  quantum  level  as  the 
characteristic  length.  This  region  is  closely  connected  into  a 
continuous  flow  region.  In  this  region,  the  shock  wave  layer 
and  boundary  layer  are  completely  mixed  together  and  we  call  it 
the  merged  layer  region.  After  the  continuous  flow  region  is 
there,  it  is  a  weak'  interaction  region  and  here  the  viscous 
region  near  the  wall  surface  and  the  non-viscous  region  of  the 
outer  surface  can  clearly  be  divided  into  two  layers.  There¬ 
fore,  the  classical  boundary  layer  theory  can  be  used  in  the 
weak  interaction  region.  Following  the  enlargement  of  the 
incident  flow's  Mach  number  M^  and  the  decrease  of  the  Reynolds 
number,  the  merged  layer  region  enlarges. 

1/2  3  1/2 

If  we  introduce  interaction  parameter  X=c  '  /(Re^  L) 

(c  is  the  Chapman  viscous  constant),  the  0  is  in  the  weak 

interaction  region  as  well  as  in  the  strong  interaction  region 
x«  &  !• 


It  has  already  been  pointed  out  that  the  shock  wave  layer  is 
completely  mixed  together  with  the  boundary  layer  in  the  merged 
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layer  region  and  therefore  use  of  the  general  boundary  layer 
theory  cannot  be  described.  In  order  to  obtain  an  even  better 
theoretical  calculation  method,  S.G.  Rubin  et  al  carried  out 
detailed  quantitative  analysis  of  steady-state  Navier-Stokes 
equations  in  different  flow  regions  and  gave  simplified  Navier- 
Stokes  equations  suitable  for  the  entire  continuous  flow 
region.  Moreover,  they  used  the  difference  method  to  numeric¬ 
ally  solve  this  set  of  simplified  steady-state  partially 
differential  equations.  The  computation  results  are  given  in 
Reference  [2] .  The  calculation  results  of  Rubin  et  al  are 
greatly  improved  as  compared  to  the  results  from  the  strong 
interaction  theory  yet  in  order  to  solve  the  set  of  simplified 
steady-state  partially  differential  equations,  it  is  necessary 
to  give  the  initial  conditions  (i.e.  the  incident  flow  boundary 
conditions)  on  a  certain  section  of  the  leading  edge  of  the  flat 
plate  and  the  corresponding  boundary  conditions.  In  this  way, 
the  computation  results  more  strongly  rely  on  the  conditions  on 
the  upstream  initial  section  and  the  differences  of  the  results 
of  different  upstream  conditions  are  very  large.  This  paper 
employs  the  method  in  Reference  [1]  and  gives  the  results  of  the 
supersonic  and  hypersonic  flows  around  flat  plates.  The  com¬ 
putation  field  is  shown  in  Fig.  2.  The  difference  schemes  is 
established  based  on  the  non-steady-state  Navier-Stokes  equa¬ 
tions  and  uses  the  initial  value  the  moment  the  uniform  incident 
flow  is  t=0.  The  upstream  boundary  is  selected  at  a  certain 
distance  on  the  leading  edge  of  the  flat  plate  and  we  also 
select  the  uniform  incident  flow  as  the  boundary  condition.  In 
this  way,  we  avoid  the  drawback  of  the  method  proposed  by  Rubin 
et  al. 


Fig.  1  Schematic  representation  of  the  network 


Fig.  2  Computational  field. 

Key:  (l)-(2.)  Incident  flow  conditions;  (3)  Three- 

point  extrapolation;  (4)  Wall  surface  conditions 
in  y=0  area. 

The  second  part  of  this  paper  uses  model  equations  to  give 
the  difference  scheme,  the  third  part  gives  the  two-dimensional 
Navier-Stokes  equations  and  their  difference  approximation 
formulas  and  lastly  gives  the  numerical  computation  results  of 
the  interaction  problems  for  the  leading  edge  of  the  flat  plate 

II.  Difference  Scheme 

For  the  purpose  of  simplification,  we  discuss  the  following 
model  equation 


In  the  formula,  v  and  a  are  constants.  Here,  we  use  the  second - 
order  derivative  term  in  space  to  simulate  the  viscous  term  in 
in  the  Navier-Stokes  equation  and  use  the  first  order  derivative 
in  space  to  simulate  the  convection  term.  It  is  necessary  that 
we  satisfy  the  steady-state  solution  of  equation  (1)  which  has 
certain  boundary  conditions.  We  use  the  following  difference 
scheme  of  Reference  [1]  for  equation  (1) 

.  ,  ..f  aJLi  —  2a£  +  ui+-l,+2  3  Ca«~u»'M)  “i.,—  “it1.  1  ,  . 

+  v  — - = -  — - —~a  2 Ax  J  (2) 

This  difference  equation  uses  second-order  accuracy  to  approx¬ 
imate  the  following  differential  equation 


an  ,  T.(  an  a-u  \ 

at  \  ax  ax- J 


In  the  equation 


At  r  (1,-2  3  )  «J1 

Ax  I  Ax  T  2  J 


u 

It  can  be  seen  that  when  =0,  equations  (3)  and  (1)  corres¬ 
pond  to  the  same  type  of  steady-state  equation.  It  is  very 
obvious  that  when  K  <  0,  the  formulation  of  the  problem  is  not 
suitable.  This  is  because  at  this  time,  as  regards  the  viscous 
term,  the  problem  of  reverse  heat  conduction  arises  and  changes 
the  characteristic  trend  of  the  convection  term  (the  charact¬ 
eristic  values  a  and  Ka  have  different  symbols).  Therefore,  in 
order  to  be  able  to  obtain  the  asymptotic  solution  when  time  t 
is  quite  large,  the  necessary  condition  is  K  ^  0,  that  is 

.^Lrvi=ii+4-i<1  c  s ) 

Ax  l  Ax  2  J 

Below,  we  will  discuss  the  stability  conditions.  Under 
K  0  conditions,  when  a  ^  0,  then  Ka  ,>  0  and  at  this  time  the 
dependent  area  of  difference  equation  (2)  always  includes  the 
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dependent  area  of  differential  equation  (3).  If  we  only  con¬ 
sider  the  convection  term,  difference  equation  (2)  is  always 
steady.  When  a  <  0,  the  -Ka  ^  0  (here,  it  is  still  necessary 
that  K  ^  0,  that  is  that  the  conditions  in  equation  (5)  are 
satisfied).  At  this  time,  in  order  to  cause  the  dependent  area 
of  difference  equation  (2)  to  include  the  dependent  area  of  the 
differential  equation,  it  is  necessary  that  A  t/A  x^  -1/Ka 
(see  Fig.  1) ,  that  is 


A*  l  A*  2  J 


(6) 


By  synthesizing  equations  (5)  and  (6) ,  we  can  obtain  the  stable 
conditions  of  difference  equation  (2) : 
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(7) 


When  /&=l/2,  the  above  stable  conditions  become 


(8) 


By  comparing  equation  (8)  and  Courant  condition  /aj 


A  t 


A  x  ^ 


^  1» 


we  can  see  that  when  taking  /9  =1/2,  the  difference  scheme  (2) 
will  doubly  widen  the  limitation  of  the  stable  conditions  towards 
step-length  ratio  At/  A  x. 


III.  The  Fundamental  Equations  and  Their  Difference  Approxima¬ 
tions 

1.  The  Fundamental  Equations 

In  the  cartesian  coordinate  system,  the  non-conservative  type 
tvo-dimensional  Navier-Stokes  equation  has  the  following  form 
after  nondirr.crsionalization  [1] 


When  there  is  nondimensionalization,  the  reference  parameters 
are  separately:  x=x/L,  y-y/L,  t=tUCD  /L,  P=p/p{0  ,  T^/T^  , 

P=p/Ufi72  “p-0  ,  u=u/Ufl?  ,  v=v/UCD  ,  ^  =M//lI<d  and  L  is  the  length 

of  the  flat  plate.  After  nondimensionalization,  the  state 
equation  is 

p  ~ySB  “T  m 


We  also  have 


’  Y  -  1 


c,-Yc, 


The  coefficient  of  viscosity  uses  the  Sutherland  relational 
formula 


In  the  formula 


T  +c. 


c,-  1  +4*-,  c.=  c.  =  114’  K 

T-  T- 


2.  The  Difference  Equations 

Using  the  difference  scheme  given  in  the  previous  section 
and  taking  t=n  A  t,  x^=i  A  x  and  y^=j^  y,  then  the  difference 
equation  corresponding  to  the  set  of  equations  in  (9)  have  the 
following  form 

UjlLzMld.  |  (a  )»  (.  /•  j  y  Uh’+i  zMj±L  =  C"  j  (13) 

At  2Ax  +  2Ajy  k  ' 


The  second-order  derivative  term  in  S.  .  (for  example,  -M—  a  -= — 

•3  3U  1 '  ^  CT  X  <7  X 

and  a  -f-)can  use  the  following  form  of  difference  approx- 

&  x  &  y 

imation 


Lt 


*x  A* 

_£_a-£H-  =  — -1— -Ca?*i./(»7*../*i 

ax  ay  4A*A y 


We  can  carry  out  similar  difference  approximations  for  other 
second-order  derivatives  of  function  u  as  well  as  the  second- 
order  derivatives  of  other  functions.  In  the  difference  equa¬ 
tions.  coefficients  a.  L.  .  and  a.  .  , .  are  separately 

1+1/2 ,3  i, 3+1/2  r  1 


a>±ui.ism-Y~(ai>i+a<±i>d 


3.  The  Boundary  Conditions 


The  wall  surface  conditions  in  the  y=0  area. 


The  non-slip  flow  conditions: 


T^Tw,  »  -  a  -  0 


Here,  the  lower  symbols  of  1  and  2  indicate  the  first  and  second 
points  of  the  adjacent  substance  surface.  The  slip  flow  condi¬ 
tions  : 

„  ,  aa 

^  TTT 


t-t,+y^y  tt(^) 


In  the  equations,  Ay  is  the  molecule's  mean  level  of  freedom. 


The  downstream  boundary  uses  the  three-point  extrapolation 
formula 

The  lower  symbol  I  indicates  the  last  point  of  the  downstream 
in  the  x  direction  (see  Fig.  2)  . 

IV.  Computation  Results 

Based  on  the  computation  method  proposed  above,  this  paper 
calculated  the  flow  interaction  problems  of  the  leading  edge  of 
the  flat  plate.  We  give  two  sets  of  computation  results  here. 
The  first  set  of  parameters  is  Ma?=20,  Tw/To=0.06,  To=9000°R  and 
RetJ?L=!4xlO  ;  the  second  set  of  parameters  is  Ma?=3,  and 

Rea?L=10^.  In  order  to  overcome  the  false  jumping  in  the  numer¬ 
ical  solution  when  there  are  super  shock  waves  near  the  leading 
edge  of  the  flat  plate,  in  the  calculations,  we  divided  the 
entire  flat  plate  into  two  sections  in  the  x  direction  and  used 
the  calculation  of  =3  as  an  example.  The  first  section  was 
-0.005  ^  x  ^  0.05  and  the  second  section  was  0.0333  ^  x  ^1.0. 

In  the  first  section,  we  used  the  step  length  Reynolds  number 
Re^x=5,  .4y/L=0.005;  in  the  second  section,  we  used 

^x/L=0.048,  y/L)  =0 . 027 ,  ( A  y/L)  m^n=0 . 02.  In  the  cal¬ 

culations  of  =20,  we  used  the  filter  method.  The  selec¬ 
tion  of  the  computational  field  and  boundary  conditions  was 
as  shown  in  Fig.  2.  We  will  now  briefly  present  the  computa¬ 
tion  results  as  follows: 

1.  Interaction  Computations  of  the  Leading  Edge  of  a  Flat 
Plate  With  Mo* -20 

In  the  computations,  the  slip  flow  conditions  were  used  on 
the  surface  and  the  density  on  the  surface  used  the  extrapola¬ 
tion  value  of  the  internal  adjacent  three  points.  The  results 
of  the  computations  are  given  in  Figs.  3-5. 


Fig.  3  Surface  pressure  distribution  (T  /T  =0.06,  =20 , 

Re  _  =4xl03 ) .  *  ° 

CO  L 

Key:  (1)  Computation  results  of  Reference  [2];  (2) 
Results  of  this  paper;  (3)  Test  results  of 
Reference  [2] . 


Fig.  4  Surface  frictional  resistance  coefficient  (T  /T  =0.06 
=20,  Refl?L=  4xl03).  W  0 

Key:  (1)  Computations  results  of  Reference  [1]; 

(2)  Results  of  this  paper;  (3)  Test  results  of 
Reference  [21. 
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Fig.  5  Temperature  profile  (Tw/To=0 . 06 ,  Mfl?=20,  Retr>L=4xlO  ). 

It  can  be  seen  from  the  computation  results  that  the  compu¬ 
tation  results  of  the  merged  layer  field  of  the  hypersonic  flow 
near  the  leading  edge  of  the  flat  plate  given  by  the  method  of 
this  paper  accords  well  with  the  test  values.  It  is  difficult 
to  give  these  kinds  of  results  when  using  the  strong  inter¬ 
action  theory.  Because  we  did  not  make  any  approximation 
assumptions  but  used  the  numerical  computation  method  to  dir¬ 
ectly  solve  the  complete  Navier-Stokes  equations,  we  more 
accurately  described  the  flow  characteristics  of  the  merged 
layer. 


2.  Computations  of  Flow  Around  Flat  Plate  With  M  =3  and 


In  this  set  of  computations,  we  took  pressure  pw=l/3 (4p1~p2) . 
In  the  formula,  lower  symbols  1  and  2  indicate  the  first  and 
second  points  closest  to  the  substance  surface.  Figures  6-8 
give  the  computation  results  of  this  set.  It  can  be  seen  from 
the  figures  that  when  Mach  number  M^,  is  not  very  large  in  the 
incident  flow,  the  influence  of  using  the  non-slip  flow  surface 
conditions  is  not  large  on  the  downstream  computations  and  that 
the  computation  results  coincide  with  the  results  given  by  the 
weak  interaction  theory. 


-***!*(  1) 

2) 

-  x*cs;tt*»*  4 
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Fig.  6  Surface  pressure  distribution  chart  (Tw  /T  =1, 
Ma>~ 3/  Re<Z?L=1°3)  ' 

Key:  (1)  Results  of  this  paper;  (2)  Computation 
results  of  Reference  [6];  (3)  Computation 

results  of  Reference  [7];  (4)  Computation 

results  of  Reference  [8] . 
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rig.  7  Surface  frictional  resistance  coefficient  (T./T  ^  -1 

Ee®L=1°  ’• 


Key:  (1)  Computation  results  of  this  paper;  (2)  Compu¬ 
tational  results  of  Reference  [6];  (3)  Computation 

results  of  Reference  [7];  (4)  Computation  results 
of  Reference  [8]. 
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Fig.  8  Velocity  profile  (*W/T0  ^=1/  M ^ =3,  Re4rL=10  ). 
(continued  next  page)  ' 


Fig.  8  (continued) 

Key:  (1)  Results  of  this  paper;  (2)  Similar  solution; 

(3)  Computation  results  of  Reference  [6] . 

Use  of  the  numerical  computation  method  to  solve  the  entire 
Navier-Stokes  equations  can  allow  us  to  better  understand  the 
flow  characteristics  of  interaction  near  the  leading  edge  of 
a  flat  plate. 
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A  NEW  METHOD  FOR  CALCULATING  UNSTEADY  SUPERSONIC  AERODYNAMIC 
FORCES  AND-  ITS  APPLICATION 

Wang  Xianxin 

Chinese  Aerodynamic  Research  and  Development  Center 
Abstract 

This  paper  proposes  a  new  method  for  calculating  supersonic 
unsteady  aerodynamic  forces.  It  combines  the  piston  theory  for 
unsteady  supersonic  flow  with  the  conical  flow  theory  for 
steady  flow  so  that  the  interaction  between  points  on  a  wing 
neglected  in  the  piston  theory  is  approximated  by  the  conical 
flow  theory.  This  method  was  used  in  this  paper  for  flutter 
calculations  of  29  cases  for  wings  with  10  different  types  of 
plane  forms.  The  results  were  compared  with  flutter  test  re¬ 
sults  obtained  in  wind  tunnels  and  they  were  in  satisfactory 
agreement.  It  was  especially  the  case  that  the  accuracy  of 
flutter  analysis  improved  by  this  method  was  noticeably  higher 
than  that  obtained  by  the  piston  theory.  This  shows  that  this 
method  is  reliable  and  feasible  for  supersonic  flutter  analysis. 
It  possesses  the  advantages  of  being  a  more  succinct  method  and 
having  higher  calculation  accuracy,  less  computing  time  and 
easier  programming. 

I.  Introduction 

The  calculation  of  the  supersonic  unsteady  aerodynamic 
forces  and  the  flutter  analysis  of  aircraft  parts  is  tedious, 
time-consuming  and  yet  indispensable  work.  At  present,  the  com¬ 
monly  used  numerical  calculation  methods  are:  the  influence 
coefficient  or  M  box  method  [1],  the  supersonic  kernel  function 
method  [2]  as  well  as  the  supersonic  doublet  grid  method  [3]  etc 
These  methods  use  the  ternary  linearized  potential  theory  as  the 
basis  and  although  they  provide  satisfactory  accuracy  when  used 
yet  none  consume  a  great  deal  of  computing  time.  An  even  more 


commonly  used  simple  and  convenient  method  is  the  piston  theory 
method  [4] .  The  practice  of  flutter  calculations  shows  the 
calculation  accuracy  of  the  piston  theory  is  sometimes  very 
good  and  sometimes  very  poor.  Therefore,  the  raising  of  the 
calculating  accuracy  of  the  piston  theory  to  develop  a  simple 
and  convenient  method  which  has  satisfactory  calculating  accur¬ 
acy  and  does  not  require  a  great  deal  of  computing  time  is  of 
real  significance. 

This  paper  proposes  a  new  method  for  calculating  super¬ 
sonic  unsteady  aerodynamic  forces.  This  method  is  formed  by 
combining  the  piston  theory  with  the  conical  flow  theory  for 
supersonic  steady  flow.  The  essence  of  the  theory  is:  the 
interaction  between  the  points  on  a  wing  neglected  in  the  piston 
theory  is  approximated  by  the  conical  flow  theory  or  we  use  the 
piston  theory  for  dynamic  frequency  revisions  and  extend  the 
conical  flow  theory  for  steady  flow  to  non-steady  conditions. 

In  this  paper,  we  carried  out  over  300  flutter  analysis  calcula¬ 
tions  of  29  conditions  of  10  wings  with  different  types  of  plane 
forms  (delta  wings,  uniform  chord  wings,  pointed  swept-back 
wings  and  trapezoidal  wings) ,  checked  this  method  and  compared 
it  with  the  flutter' test  results  obtained  in  wind  tunnels.  Com¬ 
parisons  showed  that  the  vast  majority  of  calculations  were  in 
agreement  with  the  wind  tunnel  test  data  especially  when  com¬ 
pared  with  the  flutter  analysis  carried  out  by  the  piston  theory 
wherein  the  accuracy  was  much  higher.  This  shows  that  this 
method  is  reliable  and  feasible  for  supersonic  flutter  analysis. 
It  has  the  advantages  of  being  a  succinct  method  and  having 
higher  accuracy,  less  computing  time  and  easier  programming . 

II.  General  Outline  of  the  Method 

The  interaction  of  the  points  on  a  wing  is  neglected  by  the 
piston  theory  and  the  simple  harmonic  aerodynamic  forces 
caused  by  the  airfoil  thickness  and  the  wing  surface's  unsteady 
movement  is 
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P,(X,  y,  <o;  --2Pa^ 


ico  -f-  V 
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In  the  formula,  is  the  vibration  frequency  on  the  wing  sur¬ 
face;  P  is  the  air  flow  density;  £1/  is  the  speed  of  sound;  V  is 
the  velocity  of  the  incident  flow;  Y*  is  the  specific  heat 
ratio;  h  is  the  airfoil  thickness;  W(x,y)  is  the  wing  surface's 
normal  displacement. 

The  conical  flow  theory  is  a  unique  method  for  solving  lin¬ 
earized  supersonic  steady  flows.  It  gives  the  aerodynamic  load 
on  the  wing  surface  as 

P. (x,  y)~-2pru  (2) 

When  there  is  a  subsonic  leading  edge,  disturbance  velocity  u 
is  (see  Fig.  1(a) 
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Fig.  1  Coordinate  system  and  division  of  wing  region. 

Key:  (a)  When  there  is  a  subsonic  leading  edge;  (b)  When 
there  is  a  supersonic  leading  edge;  (c)  Two- 
dimensional  region. 
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When  there  is  a  supersonic  leading  edge,  disturbance  velocity  u 
is  (see  Fig.  1  ( t») )  . 
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In  the  formulas  of  (4) ,  the  influence  of  the  quadratic  M  line 
is  neglected.  For  the  s  and  ^  in  the  formulas  of  (3)  and  (4) , 
see  Fig.  1  and  in  these  formulas  /S  -  A /m2-1 ,  M  is  the  Mach 
number,  m^/^/tgA”,  E*  (m)  is  the  second  complete  elliptical 


integral  using  ■'Sj  1-m2  as  the  modulus,  we  used  the  "calculating- 
geometric  mean  value  method"  for  calculations  and 


a0  =  $y/x 

By  using  the  P  (x,y.£b)  of  formula  (1)  and  the  P  (x,y)  of 
P  3 

formula  (2) ,  this  method  gives  the  supersonic  unsteady  aero¬ 
dynamic  forces  on  the  wing  surface  as 

p(x,  y*  <D  )-/>,(*,  y,  <a)P,(xt  y  )/Pt(y)  (5)9  F.N.l 

Here,  PQ(y)  is  the  process  numerical  function  and  its  value  is 
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(6) 


In  the  formula,  x^  and  are  separately  the  x  coordinate  of 
the  leading  and  trailing  edges  of  the  wing  surface.  In  reality. 


by  using  the  extended  theorem  of  mean  value,  PQ(y)  can  be  viewed 


as  the  pressure  of  a  certain  x  point  of  the  span  direction  y 
section.  When  each  point  of  this  section  has  this  pressure,  the 
produced  r-1  subtangent  is  equal  to  the  r-1  subtangent  produced 
by  the  conical  flow  pressure  distribution,  r  is  called  the  pro¬ 
cess  numerical  function.  It  reflects  the  calculation  method  for 
the  disturbance  quantity  of  each  point  on  the  wing  surface. 

Based  on  the  transonic  and  supersonic  law  of  similarity  to 


F.N.l  Because  the  conditions  assumed  by  the  conical  flow  theory 
are  also  approximately  satisfied  in  the  flutter  analysis  which 
is  primarily  of  the  low  order  type,  after  revising  the  piston 
theory,  the  induced  velocity  produced  by  unsteady  aerodynamic 
forces  still  approximately  satisfies  the  boundary  conditions. 
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suitably  select  similar  parameters  and  undergo  a  large  number 
of  numerical  tests,  the  formulas  for  the  given  r  are: 

For  the  subsonic  leading  edge,  m  ^  1.0 

r  -  2  +2&  (m-  6  (<L  lAt%  X  )4)  (7  a) 

For  the  supersonic  leading  edge,  m  ^  1.0 

r -(1  +4.8&/m)(  1  +3  6  (  1  -  b  ))  (7  b) 

Here, &  is  the  square  root  of  the  taper  ratio  and  A  is  the 
aspect  ratio. 

III.  Calculation  Examples 

This  method  is  combined  with  the  flutter  automated  analysis 
method  and  uses  FORTRAN  language  on  the  DJS-8  computer  to  draw 
up  the  practical  WSFA  flutter  analysis  program.  Table  1  gives 
the  circumstances  of  the  calculated  examples.  See  Fig.  2  for 
the  major  initial  data.  The  detailed  initial  data  can  be  found 
separately  in  References  [6-8].  These  examples  include  four 
different  types  of  plane  forms  as  well  as  the  situations  of 
fiat-plate  wings  and  variable  thickness  wings  and  they  are 
representative . 
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7o“?fc  =  a#  (5) 

15”¥fi3fS®«3  (g) 
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1.3,  3.0,  3.0 

1.5,  3.0,  2.5 

1.3,  1.64,  2.0,  2.33,  3.0 

1.3,  2.0,  3.0 

1.3,  2.0,  3.0 


6 

45*¥t£*52Jg«K  (g) 

1.3,  2.0,  3.0 

7 

60*^«*32fi8*  (9) 

!  1.3,  2.0,  3.0 

a 

ssfrt-aaeaflus*,  aaNACA65A004tt*mg«»  (ic 

)  1-3,  2.0 

9 

mft2.-&ag»A6Q‘,  aaNACA65Aoo4«569isa*(i]j 

.)  1.3,  2.0 

%  | 

10  i 

*a*NACA6SA002M«HS*  (12) 

1.3,  2.0 

Table  1  Examples  of  calculations. 

Key:  (1)  Example;  (2)  Description  of  wing;  (3) -(5) 
Flat-plate  delta  wing;  (6) -(9)  Flat-plate 
uniform  chord  swept-back  wing;  (10)  One-quarter 
chord  line  swept-back  angle  of  45°,  tapered 
swept-back  wing  of  NACA65A005  airfoil; 

(11)  One-quarter  chord  line  swept-back  angle  of 
60° ,  tapered  swept-back  wing  of  NACA65A004 
airfoil;  (12)  The  airfoil  is  the  trapezoidal 
wing  of  NACA65A002;  (13)  Flutter  calculation 
state  (M  number) . 
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Fig.  2  Some  initial  data  for  flutter  calculation 
examples. 

Key:  (1)  Note:  the  unit  of  size  is  cm. 

In  test  cases  4-7,  when  M=1.3,  wing 
surface  material  density 
P . =281. 16kg*  sec2/m4 . 
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The  calculation  results  are  compared  with  the  flutter  test 
data  obtained  in  wind  tunnels  and  they  are  separately  given  in 
Figs.  3-12.  The  figures  also  give  the  calculation  results  using 
the  piston  theory  method  (except  where  indicated,  they  were  com¬ 
puted  with  the  WSFA  program^ .  The  quantities  compared  were 
stiffness  height  parameter  H  and  the  flutter  frequency  ratio 

(the  ratio  of  calculated  value  and  test  value  G>  c  ).  H 

f  ,c  f  ,e 

is  the  important  dimensionless  quantity  for  determining  the 
flutter  and  it  is  given  as 

H — —  ±11 l^-swJZgnTi  +~n~Tn*)  •  (8) 

Here,  CO  is  the  first  order  torsional  frequency  of  the  wing; 

d 

Wq  is  the  weight  of  the  wing;  g  is  the  gravitational  accelera¬ 
tion. 

Figures  3-5  separately  give  the  flutter  calculation  results 
of  45°,  60°  and  70°  delta  wings.  The  test  values  were  separately 
taken  from  References  [6,7,8].  Among  them,  in  test  case  3,  the 
H  test  value  of  M=2.0  is  taken  from  the  test  curve  drawn  from  the 
test  values  [8]  but  they  are  not  taken  directly  from  test  values 
because  these  values  are  noticeably  high.  It  can  be  seen  from 
the  figures  that  aside  from  individual  points,  the  H  values 
given  by  this  method  do  not  have  more  than  6%  errors  as  compared 
to  the  test  values  and  the  errors  of  flutter  frequency  CO  ^  do 
not  exceed  8%.  The  piston  theory  method  also  has  relatively 
good  accuracy  for  the  flutter  analysis  of  delta  wings.  The 
H.  error  is  generally  about  15%  and  the  error  of  CO ^  is  about 
12%  yet  they  are  worse  than  this  method. 


Fig.  5  Comparison  of  calculated  and  test  results  for  test 

case  3. 

Key:  (1)  Test  results  [8];  (2)  This  method; 

(3)  Piston  theory  method. 

Figures  6-9  give  the  flutter  analysis  results  of  uniform-chord 
flat-plate  swept-back  wings  with  swept-back  angles  of  15°,  30°, 
45°  and  60°.  See  Reference  [6]  for  the  test  data.  It  can  be 
seen  from  the  figures  that  the  errors  of  the  H  and  dDf  values 
calculated  by  this  method  are  not  greater  than  10%,  the  indiv¬ 
idual  H  values  are  somewhat  higher  and  the  difference  reaches 
19%.  However,  under  M=1.3  conditions,  the  calculation  value 
error  of  of  each  test  case  tends  to  be  large,  20-30%.  The 
reason  for  this  is  possibly  because  in  the  initial  data  given 
in  Reference  [6],  the  test  model  material  when  M=1.3  is  not 
magnesium  again  but  was  changed  to  aluminum  alloy.  This  natur¬ 
ally  caused  the  calculated  flutter  frequency  to  have  relatively 
large  errors  (because  the  influence  of  the  change  of  the  model 
material  on  the  flutter  model  is  relatively  small) .  Further, 
in  test  case  7,  when  M=1.3  and  2.0,  aside  from  the  low  frequency 
flutter  solutions  which  occur  in  the  tests,  high  frequency  sol¬ 
utions  also  occur  in  the  calculations.  Moreover,  when  M=2.0, 
this  high  frequency  flutter  is  somewhat  lower  than  the  lower 
frequency  solution.  In  Fig.  9,  only  a  comparison  of  the  low 
frequency  flutter  solution  is  made.  The  piston  theory  method  is 


very  poor  for  analyzing  the  supersonic  flutter  of  swept-back 
wings.  All  four  of  these  test  cases  were  unable  to  calculate 
the  low  frequency  flutter  solution  in  the  tests  and  they  were 
only  able  to  obtain  high  frequency  solutions.  In  this  way,  the 
H  values  were  too  low,  reaching  more  than  50%  higher  than  the 
test  data.  The  calculated  values  of  the  piston  theory  method 
given  in  Fig.  8  were  taken  from  Reference  [4]  and  in  the  calcula¬ 
tions  it  analytically  used  polynominal  expressions  to  express  the 
test  flutter  model.  That  is,  after  this  type  of  processing  they 
were  still  only  able  to  obtain  high  frequency  solutions  and  the 
H  value  errors  were  also  15-30%.  It  is  very  clear  that  the 
accuracy  of  this  method  is  much  higher  than  that  of  the  piston 
theory  method. 
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Fig.  6  Comparison  of  calculated  and  test  results  for  test 
case  4. 

Key:  (1)  The  piston  theory  method  without  "low 

frequency"  flutter  solution;  (2)  Test  results 
[6];  (3)  This  method. 


Fig.  7  Comparison  of  calculated  and  test  results  for  test 
case  5. 

Key:  (1)  The  piston  theory  without  "low  frequency 

flutter  solution;  (2)  Test  results  [6]; 

(3)  This  method. 
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Fig.  8  Comparison  of  calculated  and  test  results  for  test 
case  6. 

Key:  (1)  Test  results  [6];  (2)  This  method; 

(3)  Piston  theory  method. 
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Fig.  9  Comparison  of  calculated  and  test  results  for  test 
case  7. 

Key:  (1)  Test  results  16];  (2)  This  method; 

(3)  Piston  theory  method. 


Two  airfoil  test  cases  with  variable  thicknesses  of 
NACA65A004,  taper  ratios  of  0.4  and  swept-back  wings  with  one- 
quarter  chord-line  swept-back  angles  of  45°  and  60°  are  given 
in  Figs.  10  and  11.  See  Reference  [6]  for  the  test  results. 


Fig.  11  Comparison  of  calculated  and  test  results  for  test 
case  9. 

Key:  (1)  Piston  theory  method  without  flutter 
solution;  (2)  Test  results  [6]; 

(3)  This  method. 

The  last  test  case  is  a  variable  thickness  trapezoidal  wing 
of  NACA65A002,  the  taper  ratio  is  0.2  and  the  one-quarter  chord¬ 
line  swept-back  angle  is  0°.  The  test  results  are  given  in 
Fig.  12.  The  calculated  results  of  this  method  are  very  good 
and  when  comparing  the  calculated  values  and  test  values  [6], 
the  H  value  errors  are  not  larger  than  5%  and  the  errors  of 
CO ^  are  not  larger  than  10%.  The  calculated  results  using  the 
piston  theory  method  for  this  type  of  non-swept-back  trape¬ 
zoidal  wing  are  acceptable  yet  the  H  differences  still  reach 
to  40-50%  and  the  CO  ^  errors  are  20%. 


Fig.  12  Comparison  of  calculated  and  test  results  for  test 
case  10 . 

Key:  (1)  Test  results  [6];  (2)  This  paper; 

(3)  Piston  theory  method. 

IV.  Conclusion 

By  comparing  the  calculated  and  test  data  of  the  flutter  of 
the  above  10  test  cases  in  29  types  of  states,  it  can  be  seen 
that  when  the  supersonic  unsteady  aerodynamic  forces  given  by 
this  method  were  used  for  flutter  analysis,  the  errors  were 
lower  than  10%  in  comparing  the  great  majority  of  obtained  H 
values  with  the  flutter  frequency  and  wind  tunnel  flutter  test 
results.  In  the  flutter  calculated  results  by  the  piston  theory 
method,  the  H  values  are  generally  larger  than  15%  and  often 
reach  40-50%.  This  is  especially  the  case  for  the  supersonic 
flutter  analysis  of  swept-back  wings.  When  using  the  piston 
theory  method  for  calculations,  it  was  not  possible  to  have  low 
frequency  flutter  solutions  appear  in  the  tests.  However,  the 
calculated  results  using  this  method  accorded  well  with  the 
test  values  which  highlights  the  advantage  of  this  method. 

This  method  is  very  succinct,  it  does  not  require  a  great 
deal  of  computing  time,  it  has  satisfactoyr  accuracy  and  it  can 
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be  used  for  model  design  units. 


I  would  like  to  express  my  deep  gratitude  to  comrade  Zhang 
Hanxin  for  his  guidance  and  help  and  to  comrade  Wang  Hongzhi 
for  his  valuable  views  proposed  for  formula  (7)  in  this  paper. 
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EFFECT  OF  A  SHEAR  LAYER  ON  THE  STABILITY  OF  AN  AXISYMMETRICAL 
EXTERNAL  COMPRESSION  AIR  INTAKE 

Zhang  Kunyuan,  Yu  Shaozhi  and  Peng  Chengyi 
Nanjing  Aeronautical  Institute 

Abstract 

This  paper  discusses  the  effect  of  six  shear  layers  with 
strengths  of  5-11%  and  variable  central  cone  external  compres¬ 
sion  air  intake  on  the  stability  of  air  intake  when  the  intake 
lip  is  in  various  positions  under  incoming  flow  M=1.72  condi¬ 
tions. 

Tests  prove  that  when  a  shear  layer  with  strength  less  than 
10%  enters  near  an  axisymmetrical  air  intake  this  does  not  cause 
separation  of  the  boundary  layer  from  the  inside  of  the  air 
intake  cowl;  shear  layers  with  strengths  up  to  11%  can  enter 
this  air  intake  at  any  position  in  the  area  of  the  inlet  and  not 
produce  buzz.  The  references  showed  that  shear  layers  with 
strengths  of  only  6-7%  can  bring  about  buzz  of  the  two  dimen¬ 
sional  air  intake  which  reveals  that  the  axisymmetric  air  intake 
has  relatively  strong  resistance  capabilities  against  the  buzz 
caused  by  the  shear  layer  intake. 

I.  Introduction 

The  stable  operating  range  of  a  supersonic  air  intake  under 
subcritical  conditions  is  an  important  area  of  its  performance. 
In  the  last  several  decades,  many  researchers  have  proposed 
their  own  views  on  this  research.  This  paper  studies  the  effect 
of  the  intake  of  a  shear  layer  on  the  stability  of  an  axisym¬ 
metrical  external  compression  air  intake. 

The  shear  layer  theory  of  Ferri  is  generally  accepted  and  it 
has  also  been  verified  by  a  great  number  of  experiments  [1] . 
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This  theory  points  out  that  when  the  shear  layer  formed  from 
the  intersection  of  shock  waves  enters  the  air  intake  near  the 
lip,  it  possibly  promotes  the  separation  of  the  boundary  layer 
from  the  inside  of  the  cowl  and  causes  air  intake  buzz.  This 
is  the  famous  Ferri  criterion. 

At  the  end  of  the  1960’s,  Fisher  et  al  of  England  success¬ 
fully  used  the  Ferri  criterion  in  a  two-dimensional  air  intake 
when  developing  the  "Harmony"  model  [2].  They  pointed  out  that 
a  shear  layer  with  strength  lower  than  6-7%  could  cause  buzz 
when  entering  the  air  intake  near  the  inside  of  the  lip.  The 
definition  of  this  strength  is  the  ratio  of  the  total  pressure 
difference  of  the  air  flow  on  the  two  sides  of  the  shear  layer 
and  the  total  pressure  of  the  upstream  air  flow  on  the  shock 
wave's  convergent  point.  The  report  also  proves  that  the  tend¬ 
ency  of  the  inner  wall  air  flow  separation  of  the  cowl  is  to 
become  stronger  with  the  strengthening  of  the  strength  of  the 
shear  layer  and  weaken  with  the  enlargement  of  the  distance  of 
the  shear  layer  from  the  cowl. 

There  is  no  doubt  that  the  research  of  Fisher  et  al  is  very 
valuable  for  two-dimensional  air  intakes.  However,  it  is 
unable  to  directly  resolve  the  analogous  problem  of  the  axi- 
symmetrical  air  intake.  These  two  types  of  air  intakes  have 
very  great  differences  in  the  reaction  level  of  the  shear  layer 
and  other  important  details.  This  paper  is  aimed  at  the  effect 
of  using  a  common  axisymmetric  air  intake  to  specially  research 
a  shear  layer  on  its  subcritical  stability. 

II.  The  Shear  Layer  of  a  Typical  Biconic  Air  Intake 

We  researched  a  double  cone  with  a  maximum  incoming  flow 
M  number  of  2.2.  The  cone  semi-angles  of  the  first  and  second 
cones  of  the  air  intake  central  cone  are  separately  17.5%  and 
25%.  When  the  incoming  flow  M  number  is  2.2  the  first  intake's 
oblique  wave  is  sealed,  when  M  ^  2.2  we  let  the  second  intake’s 


wave  be  sealed  and  when  the  M  number  is  between  1.8-2. 2, 
theoretically  the  inlet  wave  system  has  the  following  several 
situations  under  critical  and  subcritical  operations:  1. 
critical;  2.  slightly  subcritical;  3.  relatively  large  sub- 
critical  (Fig.  1)  . 


Fig.  1  Several  types  of  inlet  wave  systems  of  biconic  air 
intakes. 

We  will  now  calculate  how  large  the  produced  shear  layer ' s 
strength  is  based  on  the  structures  of  the  above  wave  systems. 

The  definition  of  the  shear  layer's  strength  Sv  is 

total  pressure  difference  on  both  sides  of  shear  layer 


v  total  pressure  of  incoming  flow  x  100% 


The  calculation  results  are  shown  in  Fig.  2.  Aside  from  this, 
distance  y  of  the  two  intake  wave's  shear  layers  from  the  lip  (or 
relative  distance  y  )  is  an  important  factor  for  determing  the  size 
of  the  effect  of  the  shear  layer.  It  is  the  function  of  the  M 
number  and  flow  coefficient  .  See  Fig.  3. 
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Fig.  2  Relationship  of  Sv  to  the  M  number. 

Key:  (1)  Head  waves;  (2)  Two  waves. 
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Fig.  3  Relationship  of  two  wave  shear  layer  distance  y  and 
flow  coefficient  <p  • 

It  can  be  seen  from  the  above  calcuations  that  in  the 
M=l. 8-2.1  range,  the  shear  layer  of  the  biconical  air  intake  has 
the  following  most  basic  situation  with  two  points: 

(1)  The  shear  layer  of  the  head  waves  is  relatively  strong 
and  its  strength  is  8.8-17.8%.  However,  this  relatively  strong 
shear  layer  cannot  enter  the  lip  near  the  critical  state  but 
can  only  enter  after  being  in  a  relatively  strong  subcritical 
state; 

(2)  The  shear  layer  strength  of  two  intake  waves  is  about 
4.8-8%.  It  can  enter  the  lip  beginning  in  the  subcritical,  its 
maximum  value  is  located  in  the  situation  whereby  flow  coef¬ 
ficient  <P<  0.9  and  at  this  time  distance  y  >0.9. 


When  analyzing  the  test  materials  published  abroad,  we 
noticed  that  the  weakest  shear  layer  which  was  able  to  cause  the 
axi symmetric  external  compression  air  intake  to  buzz  had  a 
strength  of  14.2%  [3].  Based  on  this  situation,  it  can  be  con¬ 
sidered  that  in  the  air  intake  which  we  analyzed,  a  shear  layer 


with  strengths  greater  than  14.2%  should  not  be  allowed  to  enter 
the  lip.  This  paper  focuses  research  on  the  relationship  of  a 
shear  layer  with  strength  lower  than  14.2%  but  much  higher  than 
the  6-7%  given  by  Fisher  and  the  stability  of  an  axisymmetrical 
air  intake. 

III.  Test  Model  and  Equipment 

We  used  a  single  cone  external  compression  air  intake  model 
with  variable  cone  angle  and  adjustable  cone  axial  position  and 
carried  out  shear  layer  function  tests  in  an  M-1.72  free  jet 
flow  wind  tunnel.  This  model  can  produce  shear  layers  with 
strengths  of  5.2%,  6.1%,  7.9%,  9.1%,  10.2%  and  11%  and  can  also 
enter  the  lip  under  various  distance  y.  The  model's  inlet 
diameter  is  <P  106mm,  the  wind  tunnel's  test  section  is  300x300mm 
and  the  test  which  used  the  model's  inlet  diameter  for  calcula- 
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tion  had  Re  ^2x10  -.  The  model  arranged  aloug  the  flow  has  four 
dynamic  pressure  sensors  and  total  and  static  measuring  holes 
to  measure  the  related  static  and  dynamic  parameters  under 
various  conditions.  See  Fig.  4  for  a  sketch  of  the  model. 


Fig.  4  Model  of  air  intake. 

Key:  (1)  Sensor;  (2)  Piezometric  tube. 


IV.  Test  Results 

(1)  After  a  shear  layer  with  strength  up  to  11%  enters  the 
lip,  the  static  pressure  fluctuation  of  the  air  intake's 
outlet  section  does  not  noticeably  change. 

We  separately  placed  six  different  cones  in  cone  wave  sealed 
design  cone  positions  and  during  testing  allowed  the  shear  layer 
with  six  different  strengths  to  go  from  not  entering  to  entering 
the  air  intake  near  the  air  intake  lip.  At  the  same  time,  we 
investigated  whether  or  not  the  operating  conditions  of  the 
entire  air  intake  had  any  noticeable  changes,  especially  changes 
in  the  static  pressure  fluctuation  of  the  outlet  section. 

Results  showed  that  when  these  six  cones  with  different  strengths 
entered  the  air  intake  near  the  lip,  the  air  intake  did  not  have 
any  signs  of  buzz.  The  buzz  began  in  the  air  intake  only  when 
the  positive  shock  waves  were  pushed  very  far  from  the  lip,  that 
is,  the  positive  shock  waves  began  to  have  noticeable  vibration 
and  the  outlet's  static  pressure  fluctuation  had  noticeable  en¬ 
largement.  In  brief,  a  maximum  stable  subcritical  condition 
existed.  See  Table  1. 


Table  1  (continued) 

Key:  (1)  Contents;  (2)  Cone  angle;  (3)  Shear  layer 
strength  during  testing;  (4)  Maximum  stable 
distance  (from  lip)  of  positive  shock  wave; 

(5)  Flow  coefficient  under  this  condition; 

(6)  Shear  layer  distance  under  this  condition; 

(7)  Outlet  pressure  pulsation  before  shear  layer 
enters;  (8)  Outlet  pressure  fluctuation  after 
shear  layer  enters. 


Figure  5  shows  the  curves  of  the  outlet  static  pressure 
fluctuation  of  a  24%  cone  model.  When  1  >  0.93,  although 

a  shear  layer  with  strength  of  10.7%  has  already  entered  the 
lip,  yet  the  situation  of  the  pressure  fluctuation  and  sub- 
critical  operations  do  not  have  noticeable  changes. 
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Fig.  5  Changes  in  the  pressure  fluctuation  at  outlet  section 
of  24°  cone. 

Key:  (1)  Supercritical  operation  before  shear  layer 
enters;  (2)  Stable  subcritical  operation  after 
shear  layer  enters;  (3)  Buzz. 


In  brief,  when  the  air  intake  begins  to  have  noticeable  un¬ 
steady  flow,  with  good  reason  we  consider  that  this  unsteady 
flow  was  not  created  by  the  shear  layer  but  possibly  from  the 
separation  of  the  boundary  layer  from  the  central  body. 

i 

|  (2)  The  Strongest  Shear  Layer  Does  Not  Have  Buzz  When 

Entering  the  Inlet  at  Different  Dist r^s  From  the 
l  Lip 


We  employed  two  methods  to  control  entering  distance  y  of 


the  shear  layer.  The  first  is  to  gradually  shrink  into  the 
central  cone  to  maintain  a  critical  state.  The  second  method 
is  to  gradually  push  the  positive  shock  waves  out  towards  the 
cone  tip  under  a  certain  fixed  cone  position.  In  this  way, 
distance  y  of  the  shear  layer  from  the  lip  gradually  increases 
and  gradually  approaches  the  central  cone.  We  will  now  list  a 
set  of  test  results  of  the  strongest  shear  layer  in  the 
following  table: 
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Table  2  Test  results  of  the  strongest  shear  layer  entering  the 
air  intake  at  different  distances  from  the  lip. 

Key:  (1)  Test  method;  (2)  Shear  layer  distance; 

(3)  Outlet's  static  pressure  fluctuation;  (4)  Air 
intake's  operating  condition;  (5)  Gradually 
shrinks  into  the  central  cone  to  maintain  a 
critical  state;  (6)  The  central  cone  is  fixed  and 
gradually  pushes  outwards  towards  the  positive 
shock  wave;  (7)  No  signs  of  buzz. 


This  shear  layer  did  not  produce  buzz  when  entering  our  test 
air  intake  on  various  positions  of  the  lip  section. 

We  did  similar  tests  on  shear  layers  with  five  types  of 
strength  and  the  results  were  similar.  Therefore,  no  matter 
from  what  distance  the  studied  shear  layers  with  six  types  of 
strengths  entered  the  lip,  they  were  unable  to  induce  buzz  in 
the  air  intake. 


(3)  Characteristics  of  Air  Flow  Separation  Inside  the  Cowl 


It  is  generally  considered  that  the  buzz  caused  by  the 
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entering  of  the  shear  layer  was  due  to  the  separation  of  the 
boundary  layer  inside  the  intake  cowl.  Therefore,  the  exist¬ 
ence  of  a  separation  region  brings  about  the  necessary  condi¬ 
tions  for  this  type  of  unsteady  flow.  In  this  test,  we  paid 
attention  to  the  fact  that  the  shear  layer  had  no  capability  to 
cause  the  air  intake  to  enter  buzz  operations  and  it  was  thus 
necessary  to  see  whether  or  not  these  shear  layers  had  caused 
separation  inside  the  cowl.  For  this  reason,  we  conscientiously 
researched  the  characteristics  of  the  air  flow  separation  inside 
the  cowl  after  different  strength  shear  layers  entered  the  lip 
at  different  distances.  Figure  6  gives  the  test  results.  There 
are  a  total  of  16  test  points  on  the  figure  and  when  the  strength 
was  smaller  than  10%,  7  test  points  did  not  have  air  flow  separa¬ 
tion  in  the  lip,  the  corresponding  static  pressure  fluctuation  of 
the  outlet  section  did  not  show  any  abnormal  areas  and  the  pres¬ 
sure  fluctuation  amplitude  was  only  1-2%  incoming  flow  total 

% 

pressure. 
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Fig.  6 


The  relationship  of  shear  layer  strength  S  and  distance 
y  and  the  lip  separation. 


Key:  (1)  Shear  layer  strength;  (2)  Pour  in  ink; 

(3)  Without  separation;  (4)  With  separation; 
(5)  With  serious  separation. 


When  the  strength  was  greater  than  10%,  among  9  test  points, 
lip  separation  occurred  4  times  and  among  these  the  lip  separa¬ 


tion  of  one  point  was  especially  serious.  Its  S  =10.6%,  y  &  7.4mm 
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and  yet  it  had  an  extremely  large  separation  region  and  the 
model's  outlet  static  pressure  was  still  very  small  (see  Fig.  7) 
This  shows  that  when  the  air  flow  separation  has  not  yet 
created  serious  blockage,  it  is  not  sufficient  to  induce  buzz. 


10.  ««,>■«  7.4m*  (2)24*»iia»tjspa*»;si 

Fig.  7  Outlet  pressure  fluctuations  of  24°  cone  model  without 
and  with  buzz. 

Key:  (1)  24°  cone,  S  =10.6%,  y=7.4mm,  the  outlet 

pressure  fluctuations  when  the  lip  has  serious 
separation;  (2)  Outlet  pressure  fluctuations 
when  24°  cone  has  buzz. 

Secondly,  we  noticed  that  the  shear  layer  had  the  most 
noticeable  operating  distance.  Outside  the  y  £v£  5-9mm  range, 
shear  layers  with  strengths  of  11%  are  still  unable  to  cause 
boundary  layer  separation.  In  the  tests,  the  most  noticeable 
operating  distance  y=  ^=10-18%. 

These  tests  proved  that  shear  layers  with  strength  Sv=ll% 
cannot  damage  the  axisymmetric  air  intake.  It  can  therefore  be 
anticipated  that  we  must  cause  the  shear  layer  strength  of  the 
Ferri  type  unsteady  flow  to  be  higher  than  11%  in  the  axisym¬ 
metric  air  intake. 

Naturally,  it  is  not  an  easy  question  to  resolve  how  large 
a  shear  layer  strength  a  specific  air  intake  can  sustain  and  it 
is  related  to  other  factors.  Therefore,  the  results  obtained 
from  these  tests  naturally  cannot  be  viewed  as  an  absolute 
criterion.  However,  based  on  our  test  results  and  the  published 
two-dimension  air  intake  test  results,  the  ability  of  axysym- 
metric  air  intakes  to  endure  shear  layers  must  be  stronger  than 


two-dimensional  air  intakes. 

V.  Conclusion 

(1)  The  strength  of  a  shear  layer  must  be  somewhat  larger 
than  10%  to  be  able  to  cause  separation  of  the  boundary  layer 
inside  the  cowl  of  an  air  intake.  However,  if  the  separation 
of  the  boundary  layer  does  not  produce  serious  blockage,  it  is 
not  sufficient  to  induce  buzz. 

(2)  Only  if  a  shear  layer  enters  at  a  place  a  suitable 
distance  from  the  lip  can  it  cause  the  most  marked  air  flow  sep¬ 
aration  effect.  In  these  tests,  this  distance  was  equal  to 
10-18%  of  the  inlet's  radius,  that  is,  5-9mm. 

(3)  Shear  layers  with  strengths  less  than  11%  which  enter 
axisymmetric  air  intakes  are  unable  to  cause  the  air  intake  to 
buzz.  The  axisymmetric  air  intake  can  sustain  even  higher 
strength  shear  layers  than  the  two-dimensional  air  intake. 
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COMMENT  ON  "DETERMINATION  OF  AERODYNAMIC  COEFFICIENTS  FOR  A 
REENTRY  BODY  BY  MEANS  OF  THE  KALMAN  FILTER  METHOD" 

Cai  Jinshi 

Chinese  Aerodynamic  Research  and  Development  Center 
Abstract 


The  comment  emphatically  indicates  that  a  mathematical  model  for  aero¬ 
dynamic  parameter  identification  of  a  reentry  body  bases  only  on  the  data 
of  angular  rates  and  accelerations  onboard  but  not  on  trajectory  observation  l 

data,  as  stated  in  the  paper,  it  can  only  determine  the  ratio  between  aerody¬ 
namic  coefficients  of  the  reentry  body  but  can’t  determine  the  aerodynamic 
coefficients  directly. 


In  the  latest  issue  of  the  Journal  of  Aeronautics,  Vol.  3, 

No.  3,  comrade  Jiang  Quanwei  et  al  proposed  that  in  situations 
without  orbitally  observed  data,  by  only  relying  on  data  of 
accelerations  and  angular  rates  measured  on  the  aircraft,  we  can 
determine  the  aerodynamic  coefficients  for  a  reentry  body  by 
means  of  the  Kalman  filter  method.  For  this  topic,  I  proposed 
different  views  and  discussed  them  with  comrade  Jiang  Quanwei 
and  others.  I  consider  that:  the  mathematical  model  established 
by  this  paper  to  determine  aerodynamic  coefficients  cannot  be 
identified;  when  there  is  no  orbitally  observed  data  and  we  only 
rely  on  data  measured  by  accelerations  and  angular  rates,  we  can 
only  determine  the  ratio  of  aerodynamic  coefficients  for  a 
reentry  body  but  we  are  unable  to  accurately  determine  aero¬ 
dynamic  coefficients. 

The  system  identification  theory  tells  us  that  for  the  math¬ 
ematical  model  of  a  system  identification,  if  there  is  one-to-one 
correspondence  of  the  images  from  the  parametric  space  to  the 
model's  input-output  space,  this  type  of  model  is  identifiable. 
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If  there  are  many  different  sets  of  parameters  corresponding  to 
the  same  input-output  relationship  for  the  mathematical  model 
of  certain  system  identification,  then  this  model  is  not 
identifiable. 

The  parametric  space  of  this  paper's  system  identification 
mathematical  model  is  composed  of  the  unestimated  value  of 
aerodynamic  coefficient  C  and  the  state  variable  initial  value 
(—o,  Hq) .  Rudder  angle  £  is  known  and  is  the  system  input. 

Because  this  model  relies  on  accelerations  and  angular  rates  to 

A/  /V 

carry  out  identification,  the  output  of  the  system  is  (N  ,  N  , 

x  y 

N  ,  P»  <3  and  r)  •  This  is  easy  to  verify. 

For  this  system,  numerical  group  A(C,  X^and  Hq)  and  data 
B(Cexp(  A  I^/k^X^,  Hq+  4Hq)  of  the  parametric  space  correspond 
to  a  similar  input-output  relationship  and  therefore  the  group 
of  parameters  is  not  identifiable. 

By  only  relying  on  data  of  acceleration  and  angular  rates, 
we  cannot  determine  that  the  aerodynamic  coefficients  have 
physical  significance.  Aircraft  movement  is  determined  by  the 
force  and  moment  acting  on  the  aircraft;  aircraft  movement  para¬ 
meters  (accelerations  and  angular  rates)  can  also  in  turn  deter¬ 
mine  the  force  and  moment  on  an  aircraft.  However,  aerodynamic 
force  and  moment  are  a  function  of  the  product  of  the  aerodynamic  coefficient 
and  the  dynamic  pressure  Q.  For  aircraft  with 

similar  initial  states,  different  aerodynamic  coefficients  only 
require  that  it  have  the  same  product  as  the  dynamic  pressure  to 
be  able  to  have  similar  aerodynamic  force  and  moment  and  thus 
have  the  same  accelerations  and  angular  rates.  Therefore,  by 
only  relying  on  accelerations  and  angular  rates,  if  there  is  no 
new  additional  information  to  separate  out  the  dynamic  pressure, 
then  there  will  be  no  way  to  determine  the  aerodynamic  coef¬ 
ficient.  This  is  the  physical  nature  of  things  and  cannot  be 
resolved  by  a  mathematical  method  (e.g.  Kalman  filter). 
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By  only  relying  on  the  data  of  accelerations  and  angular 
rates,  we  are  only  able  to  determine  the  ratio  of  the  aircraft's 
aerodynamic  coefficients.  Naturally,  if  we  are  able  to  deter¬ 
mine  a  certain  aerodynamic  coefficient,  we  can  also  determine 
other  aerodynamic  coefficients. 

From  formula  (9)  of  the  paper  we  can  obtain 

Q=‘CN.,  +  Xlm(g1+rt)'}/(c,as)s(Nxg)/(c,Bs') 

By  substituting  the  above  formula  into  formulas  (8)  and  (9) 
of  the  paper,  we  can  obtain  a  new  set  of  state  equations  and  a 
set  of  observational  equations.  Its  undetermined  aerodynamic 
coefficient  is 


The  state  equations  and  observational  equations  of  this 
mathematical  model  do  not  include  dynamic  pressure  Q  and  height 
H.  Therefore,  it  is  not  necessary  to  have  orbital  parameters 
and  by  using  the  maximum  likelihood  method  or  Kalman  filter 
method  we  can  identify  the  ratio  of  the  aerodynamic  coefficients 
but  cannot  identify  the  aerodynamic  coefficients. 

This  paper  used  the  generalized  Kalman  filter  method  and  used 
the  aerodynamic  coefficient's  filter  convergence  value  as  the 
determined  aerodynamic  coefficient  value.  Because  of  the 
unidentif iability  of  these  coefficients,  the  initial  filter  value 
of  the  filter  convergence  value  and  state  parameters  as  well  as 
the  initial  covariance  value  are  related.  Simulated  calculations 
show  that  by  changing  these  initial  values,  the  determined  aero¬ 
dynamic  coefficients  will  also  be  different.  It  is  only  neces¬ 
sary  for  these  initial  values  to  slightly  allow  modification  and 
then  the  errors  of  the  main  determined  aerodynamic  coefficients 


can  reach  above  20%.  Therefore,  when  resolving  in  situations 
without  orbitally  observed  data  spoken  of  in  this  paper,  it  is 
not  precise  to  only  rely  on  the  data  of  accelerations  and 
angular  rates  to  determine  coefficients. 


THIRD  NATIONAL  CONFERENCE  OF  THE  AERONAUTICAL  SOCIETY  OF  CHINA 


The  Aeronautical  Society  of  China  held  its  Third  National 
Conference  in  Xian  from  March  24  to  28  in  1983.  Nearly  300 
representatives  from  aeronautical  and  aerospace  production, 
utilization,  scientific  research  and  mathematics  departments 
from  various  districts  of  China  attended  the  conference.  The 
conference  listened  to  and  deliberated  on  the  work  report  of  the 
Second  Council,  discussed  the  passage  of  the  new  constitution 
and  employed  the  secret  ballot  method  to  select  a  new  council. 
Wei  Wenmei  was  elected  as  director  and  Gao  Zhenning,  You  Jiang, 
Yao  Jun,  Wang  Wanlin,  Li  Ming,  Zhuang  Fenggan,  Cao  Chuandiao 
and  Zhang  Azhou  were  elected  as  deputy  directors  and 
Wang  Nanshou  was  elected  as  secretary-general. 

The  standing  council  approved  inviting  Cao  Lihuai,  Mei 
Jiasheng  and  Shen  Yuan  as  advisers  to  the  Aeronautical  Society 
of  China. 


AERONAUTICAL  SOCIETY  OF  CHINA  HOLDS  ITS  "SECOND  ACADEMIC 
EXCHANGE  CONFERENCE  ON  ENGINE  STRUCTURE,  STRENGTHS  AND  VIBRATIONS 

The  Aeronautical  Society  of  China  convened  its  "Second 
Academic  Exchange  Conference  on  Engine  Structure,  Strengths  and 
Vibrations"  in  Wuxi  in  Jiangsu  Province  from  April  21  to  25  of 
1983.  The  conference  representatives  were  over  140  persons  from 
a  total  of  46  units  of  the  airforce,  navy,  civil  aviation,  the 
Ministry  of  the  Aerospace  Industry,  the  Ministry  of  the  Aviation 
Industry,  the  Ministry  of  the  Machinery  Industry,  factories 
attached  to  the  Chinese  Academy  of  Sciences  system  and  related 
institutions  of  high  education  from  throughout  China.  The  con¬ 
ference  received  106  academic  papers  and  reports,  73  of  which 
were  read  at  the  conference  and  there  were  also  two  specially 
invited  reports. 

The  papers  and  reports  read  at  the  conference  were  all  con¬ 
cerned  with  the  newest  achievements  in  application  research  in 
the  fields  of  engine  structural  design,  strength  and  vibration, 
materials  and  material  strength  etc.  They  prominently  showed 
the  advancements  in  recent  years  in  the  fields  of  aviation 
engine  structural  integration  and  aviation  engine  developed 
programs,  engine  load  spectrum  and  lifetime  calculation  research, 
engine  materials  and  materials  research  etc.  Related  experts 
were  also  invited  to  the  conference  to  present  comprehensive 
reports  on  the  state  of  present  domestic  and  foreign  technical 
advancements  and  these  experts  were  welcomed  by  the  represent¬ 
atives. 

During  the  conference,  separate  symposiums  by  factories, 
research  institutes  and  institutions  of  higher  learning  were 
held  on  the  problems  surrounding  future  developments  in  engine 
structures,  strength  and  vibrations  and  a  special  symposium  on 
high  temperature  fatigue  strength  was  also  held.  During  the 
conference,  the  formal  establishment  of  "Engine  Structure, 
Strength  and  Vibration  Special  Group  of  the  Special  Committee 
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on  Dynamics  of  the  Aeronautical  Society  of  China  was  declared 
and  the  first  working  conference  of  the  group  was  convened. 
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THE  AERONAUTICAL  SOCIETY  OF  CHINA  HOLDS  ITS  "ACADEMIC  CONFERENCE 
OF  AVIATION  METAL  CEMENTING" 

The  Aeronautical  Society  of  China  convened  the  "Academic 
Conference  of  Aviation  Metal  Cementing"  in  Tongan  County  in  Fujian 
Province  from  April  11  to  15  of  1983.  The  conference  represent¬ 
atives  included  over  100  persons  from  47  units  of  the  airforce, 
navy,  civil  aviation,  the  Ministry  of  the  Aerospace  Industry, 
the  Ministry  of  the  Aviation  Industry,  the  Ministry  of  Petrifica¬ 
tion,  factories  attached  to  the  Chinese  Academy  of  Sciences 
system,  research  institutes  and  related  institutions  of  higher 
education.  The  conference  received  a  total  of  48  papers  and 
reports  and  among  these  10  were  read  at  the  plenary  session. 

The  papers  and  reports  read  at  the  conference  were  all  con¬ 
cerned  with  the  newest  achievements  in  application  research 
and  experimental  research  in  the  fields  of  cementing  structural 
design,  cementing  techniques,  adhesive  materials  and  detection 
techniques  etc.  They  basically  reflected  the  present  advanced 
level  of  China's  metal  cementing  technology.  Related  experts 
were  also  specially  invited  to  the  conference  to  present  com¬ 
prehensive  reports  on  the  state  of  present  domestic  and  foreign 
technical  advancements  in  cementing  technology.  They  also 
showed  scientific  and  technical  films  related  to  technical 
advancements  in  cementing  technology  and  thus  invigorated 
scientific  thoughts,  widened  scientific  outlooks  and  this  was 
inspiring  to  the  representatives  and  received  their  unanimous 
favorable  comments.  The  conference  also  carried  out  specialized 
discussions  of  present  problems  existing  in  the  development  of 
China's  metal  cementing  technology  and  future  development 
trends  and  they  also  proposed  consultative  recommendations  to 
related  departments. 


